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We present an overview of recent developments in numerical relativity studies of higher 
dimensional spacetimes with a focus on time evolutions of black-hole systems. After a 
brief review of the numerical techniques employed for these studies, we summarize results 
grouped into the following three areas: (i) Numerical studies of fundamental properties of 
, black holes, (ii) Applications of black-hole collisions to the modeling of Trans-Planckian 

i scattering, (iii) Numerical studies of asymptotically anti-de Sitter spacetimes in the con- 

text of the gauge-gravity duality. 
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1. Introduction 

Numerical modeling of dynamic spacetimes containing black holes (BHs) in the 
framework of Einstein's theory has been motivated for a large part of its nearly 50 
year long history by the calculation of the gravitational wave (GW) signals generated 
by BH binary systems. These efforts culminated in 2005 in the first evolutions of 
BH binaries through inspiral, merger and ringdownP^ Over the following years, 
numerical relativity (NR) simulations have been used in the construction of GW 
template banks 4 " 8 fundamental for signal searches in GW detectors such as LIGO, 
VIRGO and GEO60GPE 1 and have also been studied in the context of GW data 
analysisP^^ The Ninjifl^ and NRAR^ collaborations in particular combine the 
resources of various NR groups with the analytic relativity and gravitational-wave 
data analysis community. NR studies have furthermore provided valuable insight 
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into a variety of astrophysical scenarios involving BHs as for example the study of 
gravitational recoilpH^ID j e t phenomena^^ and electromagnetic counterparts for 
multi-messenger astrophysics p3M] p or reviews of these subjects we refer to Refs.[29l 
[301 EH [321 [331 

In recent years, however, NR has also emerged as a valuable tool to explore high- 
energy physics scenarios.^ The purpose of this article is to provide an overview 
of the most recent developments in these areas which, for reasons of clarity, we 
divide into three categories: (i) numerical studies of fundamental properties of BH 
spacetimes, (ii) applications to the modeling of Trans-Planckian scattering and (iii) 
the gauge gravity duality. This division is not entirely rigorous and we will encounter 
scenarios of relevance for more than one of these major topics. For example the 
stability of rotating BHs in higher dimensions is discussed in Sec. [3] on fundamental 
properties of BHs but likely also affects the signature in the conjectured formation of 
mini BHs in Trans-Planckian scattering events. Wherever present, we will emphasize 
such significance across our somewhat arbitrarily chosen distinction. 

NR studies in the context of high-energy physics and fundamental properties of 
classical general relativity frequently involve BHs in higher dimensional spacetimes. 
In some cases, however, valuable insight can be inferred as well from the more 
traditional case of spacetimes with "3+1" dimensions (three spatial and one time- 
like). Somewhat in breach of its title, this work is therefore not entirely limited to 
the numerical investigation of higher-dimensional BHs. As complementary reading 
to this article, we note review articles including NR applications in high-energy 
physics by PretoriusJ^ Cardoso et al¥^ and Yoshino & ShibataP^^ An extended 
discussion of BHs in higher dimensions is given in [351 

The remainder of this article is organized as follows. We begin with a brief 
summary of the NR techniques in Sec. [5] The above mentioned main applications 
of NR to higher-dimensional BHs will be reviewed in Sees. G2 [H and [SJ In each of 
these we will provide a brief motivation of the corresponding studies and list more 
specialized review articles for further reading. We will conclude in Sec. [6] Unless 
noted otherwise, we use units such that the speed of light c = 1. 

2. Numerical techniques in higher dimensions 

Formulations of the Einstein equations suitable for NR can be classified into two 
main groups. The characteristic Bondi-Sachs^"^ formalism is based on the char- 
acteristic surfaces of the vacuum Einstein equations. This formalism leads to a 
natural hierarchy of the equations which reflects the isolation of the gravitational 
degrees of freedom; in the four-dimensional case, for example, one obtains one trivial 
equation, three complementary equations, four hypersurface equations and two evo- 
lution equations for the degrees of freedom. A difficulty arising in the characteristic 
approach is the potential formation of caustics and the ensuing breakdown of the co- 
ordinate system. For the classical general-relativistic two-body problem, the inspiral 
and merger of two BHs, a robust solution for this problem has not yet been found. 
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Numerical relativity in higher dimensions 3 

Nevertheless, characteristic methods have been employed with great success in var- 
ious codes studying configurations with additional symmetries, most notably in the 
series of studies by Chesler & Yaffe^3EH on thermalization of quark-gluon plasmas 
through the Anti-de Sitter /Conformal Field Theory (AdS/CFT) correspondence. 
Characteristic evolutions are reviewed in detail by WinicourPH 

The majority of work on the inspiral and coalescence of compact binaries, in- 
cluding the 2005 breakthroughs, however, has been performed in 3+1 approaches. 
Here spacetime is decomposed into a one-parameter family of three or, more gen- 
erally, d — 1 dimensional, spatial hypersurfaces. The canonical 3+1 split has been 
developed by Arnowitt, Deser and Misner^ (ADM), reformulated by YorlP^ anc j 
leads to a system of six second-order evolution equations for the induced spatial 
metric and four constraint equations, the Hamiltonian and momentum constraints. 
Introduction of the extrinsic curvature as an evolution variable converts the evolu- 
tion equations into a first-order-in-time system. By performing a conformal rescaling 
of the spatial metric, a decomposition of the extrinsic curvature into trace and a 
trace-free part and evolving the contracted Christoffel symbol as separate fields, 
Baumgarte, Shapiro, Shibata and Nakamura^^l have developed a strongly hyper- 
bolic formulation^ commonly referred to as "BSSN" which has been used with 
great success by many numerical groups. The ADM/BSSN approach to numerical 
relativity is discussed in great detail in [50l EU E21 

An alternative Cauchy-type formulation has been derived from the Einstein 
equations in harmonic gauge where they take on a manifestly hyperbolic form. 
The harmonic gauge can be generalized to arbitrary gauge choices by introducing 
source functions for the coordinate conditions .ESUHl More details of this generalized 
harmonic gauge (GHG) formulation can be found in[5SJ[SSl[S71 The GHG and BSSN 
methods have been the foundation of the numerical relativity breakthroughs of 2005 
by Preterms^ and the Brownsville and Goddard groupsP^ 

NR applications in higher-dimensional spacetimes have so far been restricted 
to configurations with symmetry assumptions. This is largely a consequence of the 
vast increase in computational resources required by each additional dimension. 
A reduction of higher-dimensional general relativity to an effective 3+1 or lower- 
dimensional computational domain has been achieved in various ways. The most 
evident approach is to directly impose the symmetry together with the coordi- 
nate conditions on the line element. For example, Sorkin & Orerl^l describe the 
(i-dimensional, spherically symmetric spacetime in double-null coordinates (it, v) by 
the line element 

ds 2 = —a(u, v) 2 du dv + r(u, v) 2 dri 2 d _ 2 > (1) 

where d^ld-2 is the metric on a d— 2 dimensional unit hyper sphere and r is the areal 
radius. The time evolution is then determined by the Einstein equations derived 
from this metric; see their Eqs. (3)- (6). 

Alternatively, a general procedure for a dimensional reduction under the assump- 
tion of one or more Killing vectors is given by Geroch.^l By using this procedure, 
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the d-dimensional Einstein equations with SO{n) symmetries can be converted to 
an effective 3+1 scheme coupled to scalar or vector fields which represent the cur- 
vature of the extra-dimensions.^^ Such a dimensional reduction is particularly 
attractive because it provides a straightforward way to apply the successful numer- 
ical techniques developed in 3+1 with relatively minor modifications to the study 
of higher dimensional spacetimes. 

A different way to achieve the same goal is based on the cartoon method^ 
which employs interpolation of grid functions and/or trading of spatial derivatives 
to evolve symmetric spacetimes on a reduced Cartesian grid. This approach has 
been modified and extended to higher-dimensional applications by Pretorius^ and 
Yoshino & Shibata 37 63 and has proven a rather robust tool to study BH collisions 
and instabilities. 

The generation of initial data in higher-dimensional spacetimes relies to a signif- 
icant extent on similar techniques developed for the 3+1 case; see Cook's^ review. 
In particular, it has been shown by Yoshino et alW^ that the existence of analytic 
solutions of Bowen-YorlpSJ type for the momentum constraints carries over to gen- 
eral relativity in d dimensions. Solving the Hamiltonian constraint for the conformal 
factor can then be achieved by numerical methods similar to those developed for 
3+1 dimensions. By using this procedure, initial data for BH binaries with boost 
have been obtained in Refs. |65j EH For specific numerical studies, analytic initial 
data is available in closed form, as for example the d dimensional generalization of 
Brill-Lindquistr^ data describing a BH binary at the moment of time symmetry. 
Analytic BH solutions can also be employed to construct approximate solutions 
to the Einstein constraints by superposition of either the solutions themselves or 
boosted variants thereof. 

For more details on the numerical methods, the construction of initial data and 
gauge choices, we refer the reader to the publications reviewed in the following 
sections. 

3. Fundamental properties of black holes 

BHs play a central role in the applications of NR to higher-dimensional spacetimes 
discussed in Sees. |4] and O Quite aside from these applications, BHs provide a unique 
probe for deepening our insight into the gravitational theory and understanding in 
which regards 3+1 dimensions represent a special case. BHs in higher dimensions 
indeed reveal a much richer phenomenology than their four-dimensional cousins; 
cf. Emparan & ReallP™ 

Asymptotically flat, stationary solutions of the Einstein equations in vacuum 
in 3+1 dimensions belong to the Kerr family and describe single BHs of spherical 
topology characterized by two parameters, the mass M and angular momentum 
parameter a.^^ The stability of these solutions has been studied extensively; see 
for example Refs. [74j [75] and Dafermos^ review. Such classical BH solutions are 
expected to be stable in the region outside the horizon, but admit instabilities in the 
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interior.^ Penrose's cosmic censorshi]P^^ furthermore conjectures that collapse of 
ordinary matter does not lead to the formation of naked singularities; instead such 
singularities are hidden from outside view through an event horizon. 

BHs in higher dimensions are solutions to the d-dimensional Einstein equation^ 

G U v = Ry.i> — -^9uuR = 8TTGT uly , (2) 

for the vacuum case T^ v = 0. 

Solutions to these equations can be generated straightforwardly from extending 
the four-dimensional Schwarzschild metric with one or more flat dimensions. These 
black strings or black branes are susceptible to the Gregory-Laflamm^ilSlI (GL) 
instability related to the vastly different length scales in the horizon geometry. These 
solutions furthermore suggest a modification of Thome's^ Hoop conjecture; if an 
amount M of mass-energy is compressed into a volume such that its circumference 
is below its Schwarzschild radius in every direction, a BH forms. The existence of 
black branes indicates that such a concentration of energy is sufficient in a subset 
of the spatial directions. 

The d-dimensional generalization of the spherically symmetric Schwarzschild 
solution is given by the Tangherlini^ metric 

ds* = - (l - df + (l - - 1 dr* + r>dtf d _ 2) (3) 

where the constant // determines the BH mass M = (d — 2)f2<2-2/V(167rG). Here, 
£ld-2 denotes the surface area of the d — 2 dimensional unit hypersphere. Note that 
mass is now of dimension length^ -3 . Schwarzschild solutions in d dimensions have 
been found to be stable against linearized gravitational perturbations ! 84 ^ 85 ! 

In contrast to their spherically symmetric counterparts, the class of rotating 
BH solutions exhibits a vastly richer structure in higher dimensions. First, rotation 
is now possible in more than one plane, giving rise to more than one rotation 
parameter. Second, the radial dependence of the Newtonian potential is altered 
to l/r d ~ 3 whereas the centrifugal potential maintains its 1/r 2 character. This d- 
dependent character of the Newtonian potential manifests itself in the absence of 
stable circular geodesies for d > 4.HSIH3 

Quite remarkably, Myers & Perry^ found an exact solution describing BHs 
rotating in all possible planes. In d = 4 and d = 5, BHs rotating in one plane with 
rotation parameter a > M and a 2 > /i respectively represent naked singularities. 
In d > 6, in contrast, there exists a horizon for all values of a. Furthermore, in 
the limit of infinite rotation rate, this horizon becomes increasingly flattened along 
the rotation plane^ and is therefore expected to be GL unstable. This has been 
confirmed in linearized studies by Dias et alW^H In the case of rotations in more 
than one plane, the existence of horizons in ultra-spinning BH spacetimes requires 



a see Sec. 3.1 in 1701 for a discussion of Newton's constant G 
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at least two (one) spin parameter to vanish for an odd (even) number of dimensions 
d. 

Rotation in higher dimensions also leads to BH solutions of non-spherical topol- 
ogy. Black rings are supported against contraction by the centrifugal potential; 
analytic solutions in d = 5 have been found in [92l [93] Thin rings behave locally like 
boosted black strings^ which are GL unstable 9 -^ while fat rings are expected to be 
unstable to variations of their radius! 9 ^ Solutions for black saturns have been con- 
structed in 1971 and demonstrate the existence of stationary vacuum multiple-BHs in 
higher dimensions. 

The generalization of the AdS-Schwarzschild solution^ is given by 

^ - - (l - £s + i ) + (l - £s + i) 1 ** + rW d _ 2 . (4) 

In the limit fj, — > 0, this solution reduces to the AdS spacetime with curvature ra- 
dius L in global coordinates; cf. Sec. [5] Its mass relative to the empty AdS solution 
is given b y 99 '' 100 ' M = (d — 2)f^_2/J (16ttG). This solution is stable against lin- 
earized gravitational perturbations!^' The four-dimensional Kerr- AdS spacetime 
(see e. g. Eq. (2.1) in I102p has been found by Cartei^^ and generalized to higher 
dimensions in Refs. 11041 11051 11061 Kerr- AdS has been found to be unstable against 
linear perturbations if the BHs are small ES3tll2l 

One of the first and most influential NR results has been obtained in Chop- 
tuik'iJiiil seminal study of the collapse of spherically symmetric massless scalar 
fields minimally coupled to Einstein gravity in four-dimensional, asymptotically flat 
spacetimes. By evolving various one-parameter families of scalar pulses, he identified 
critical behaviour as the parameter p which characterizes the gravitational interac- 
tion strength of the field approaches a critical value p* . For p > p* , the field collapses 
to a BH and for p < p* it disperses to infinity. Furthermore, near-critical field con- 
figurations exhibit universal behaviour in the strong field limit: (i) BHs which form 
have a mass M ~ \p — p*| 7 with a universal constant 7 0.37. (ii) Advancing the 
evolution from a time t to t + A, the field profile is recovered up to a "zoom-in" by 
a factor e A . The numerical study reveals A to be a universal constant of about 3.4. 
For subcritical configurations, Garfinkle & Comer Duncaif^ found a similar scaling 
of the maximal scalar curvature in the spacetime i? ma x ~ \p — P*\ 21 ■ Continuously 
self-similar solutions were found by Pretorius & Choptufif^ for scalar fields in 2+1 
dimensional asymptotically AdS spacetimes with a mass scaling characterized by 
7/2 = 1.2. Sorkin & OrerP^ generalized Choptuik's result to higher dimensions by 
evolving scalar fields up to d = 11 dimensions. Their results indicate that 7 reaches 
a maximum and A a minimum around d ~ 11 . . . 13. An extended review of critical 
collapse phenomena is given by Gundlach & Martfn-GarcfaP^I 

More recently, critical-collapse studies have been extended to asymptotically 
AdS spacetimes in 3+1 and higher dimensions. In a remarkable study, Bizoh & 
Rostworowski found evidence suggesting that the 3+1 AdS is unstable to BH for- 
mation under arbitrarily small perturbations. By evolving spherically symmetric 
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scalar fields they recover Choptuik's results for large initial field amplitudes. For 
configurations below this critical value, however, the AdS boundary substantially 
modifies the outcome. In contrast to asymptotically flat spacetimes, spatial infinity 
in AdS is reached in finite time by massless fields and reflected back onto the origin. 
As the initial amplitude is reduced below a critical value eo, the scalar pulse forms a 
BH upon its second implosion on the origin. Further reduction of the amplitude leads 
to a second critical amplitude t\ and this pattern repeats itself with no indication of 
a threshold amplitude for BH formation; cf. their Fig. 1. For each critical amplitude, 
they furthermore recover Choptuik's scaling law with 7 = 0.37. They interpret this 
behaviour as a resonant mixing of modes which transfers energy from low to high 
frequencies. This study has been generalized to higher-dimensional AdS spacetimes 
by Jalmuzna et aiP^ suggesting that AdS is unstable to BH formation for generic 
spacetime dimensions. Presumably, this result has eluded a similar study by Garfin- 
kle & Pando Zayasp^ because of insufficient length of their numerical simulations 
for smaller field amplitudes. In a similar investigation using complex scalar fields, 
Buchel et alW^ reproduce the instability of AdS and the transfer of energy from low 
to high frequencies. In consequence, the width of initially weak pulses narrows in 
each reflection cycle and eventually collapses to a BH; cf. their Fig. 5. Garfinkle et 
a ^jii8j j^vg monitored the time of BH formation and the horizon radius and find the 
amplitude of the scalar field to have a stronger influence on the outcome compared 
with the width of the pulse. The same type of instability to BH formation has been 
found by MaliborskfS^l f or a Minkowski spacetime enclosed inside a reflecting wall, 
indicating that the global structure plays a major role for the effect. Perturbative 
studies support the numerically observed instability of AdS but the question of 
the generic stability properties of asymptotically AdS spacetimes, as for example 
those containing BHs or boson stars, remains under study 

The instability of black strings has been studied numerically in a sequence of 
papers by Lehner & Pretorius and collaborators LLSSHIMI see a i so |125I It had been 
known since the 1990s that black strings are subject to the GL instability, but the 
eventual fate of the string remained unclear. In Ref. 11241 Lehner & Pretorius found 
evidence supporting indications by earlier work that the string evolves to a sequence 
of spherical BHs connected by thin string segments which themselves are subject 
to the GL instability, resulting in a self-similar cascade reaching zero string width 
in finite asymptotic time. This behaviour shows striking similarity with satellite 
formation in the flow of low- viscosity fluids. The eventual bifurcation of the horizon 
resulting from the cascade implies formation of a naked singularity P^ Because no 
finetuning is required to trigger the instability, the result constitutes a violation of 
the cosmic censorship conjecture without "unnatural" assumptions about the initial 
data. In contrast to the higher-dimensional case, NR has as yet not observed any 
such violation of cosmic censorship in 3+1 dimensions. In a recent study^S BH 
collisions in asymptotically de Sitter spacetimes have been found to comply with 
censorship; BH binaries with a combined mass exceeding the inverse of the Hubble 
constant do not merge for any initial separation provided the initial data do not 
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contain a naked singularity. 

We have already mentioned the linear instability of rapidly rotating Myers- Perry 
BHs. The stability of rotating BHs has been studied numerically by Yoshino & Shi- 
bat£p^2Sl w ] 10 evo ive single holes with a bar-shaped perturbation imposed on the 
conformal factor in 5 < d < 8 dimensions. In all cases, they observe BHs spinning 
above a threshold dimensionless spin-parameter j cr it to spontaneously emit gravi- 
tational waves and settle down into a subcritical configuration. Expressed in terms 
of the ratio of polar to equatorial circumference C p /C e of the apparent horizon, the 
onset of instability occurs for nearly identical values for d > 6, namely 0.65, 0.68 
and 0.67, respectively, for d = 6, 7 and 8. In contrast, they obtain a significantly 
stronger distortion threshold 0.38 in d = 5 indicating a somewhat special status of 
d — 5; we recall in this context the absence of an upper limit on the spin in d > 6. 
BHs spinning below the critical rate, on the other hand, show no signs of instability 
in any of their simulations. The effect of this instability on the evaporation of BHs 
conjectured to be produced in trans-Planckian scattering (cf. Sec. |4|) depends on 
the timescale of the instability and, hence, on the spin parameter of the formed BH; 
for details see the discussion in Sec. VI D in Ref. 129. 

Even though this review is primarily concerned with time evolutions, we note 
that numerical techniques are also applied for finding static and stationary solutions 
to the d dimensional Einstein equations in vacuum, as for example the construction 
of four-dimensional rotating BH solutions in a cavitjff^SI or asymptotically AdS 
spacetimes in d = 5 with a conformally Schwarzschild boundary This topic is 
reviewed in detail by Wiseman L122J 

4. Trans-Planckian scattering 

The standard model of particle physics provides an exceptionally successful descrip- 
tion of subatomic particles and their interactions via the electromagnetic, weak and 
strong forces. In spite of its success, however, there remain important unanswered 
questions, as for example the unknown nature of dark matter. In the context of 
this section, the most important open issue is the large discrepancy between the 
electroweak energy scale 246 GeV and the grand unification scale ~ 10 19 GeV. 
This so-called hierarchy problem manifests itself in the extraordinary weakness of 
the gravitational interaction relative to the other fundamental forces; the weak in- 
teraction is about 32 orders of magnitude stronger than the gravitational one. At 
energies comparable to the grand unification energy, on the other hand, gravity is 
expected to become comparable in strength to the other interactions. 

A particularly intriguing solution to the hierarchy problem is the ADD 
modelp^HHH which proposes large (relative to the four-dimensional Planck scale 
Mpi) extra dimensions of radius ~ R such that the fundamental Planck scale Mpi.d 
may be as low as ~ 1 TeV. Consider for illustration the Poisson equation in d — 1 
spatial dimensions for a point source A$ = AitGdM , with the d dimensional gravi- 
tational constant Gd = h d ~ 3 c 5 ~ d Mp^ d i and the source mass M. From Gauss's law 
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we obtain in spherical coordinates 
I d := / A$dV = f V$dA = d r <S>r d - 2 n d _ 2 = ^ ^J^V ' 2 = 4nG d M . (5) 

J J 1 ( 2 J 

At large distances r ^> R, the extra dimensions are not accessible to the gravita- 
tional field and merely contribute an overall volume factor such that I d ~ R A Ii 
whereas for r < R, the potential is determined by the d — 1 dimensional Poisson 
equation. Ignoring geometrical factors of order Oil), we thus obtain 

G d M G d M 

fOT r <R, $^-_iL_ for (6) 

At large distances we thus recover the usual 1/r potential with a four-dimensional 
coupling constant G 4 = Gd/R d ~ A which implies for the energy scales Mp^ ~ 
MpfJiJ d-4 . For example, a fundamental Planck mass Mv\,d of the order of the 
electroweak scale requires R ~ 10 13 cm in d = 5 (ruled out experimentally) or 
i? < 1 mm for rf = 6. In a different version of the model proposed by Randall 
& Sundrump^E^ a finite length scale is introduced through a warp factor into 
otherwise infinite extra dimensions. 

An intriguing consequence of an effective Planck mass much below its four- 
dimensional value ~ 10 19 GeV is the possibility of BH formation in parton-parton 
collisions at the Large Hadron Collider (LHC)P^'i 3 - 9 -' se e also the reviews by 
CavagligP^Sl and KantiP^Once formed, such mini BHs are expected to evaporate in 
four stagesfHB (i) a balding phase during which the BH sheds all multipoles except 
for mass, spin and charge, (ii) a spin-down and (iii) a Schwarzschild phase during 
which the BH looses first its spin and then its mass via semi-classical Hawking 
radiation and (iv) the Planck regime as the BH mass approaches the Planck mass 
which is described by an as yet unknown theory of quantum gravity. 

Of particular interest in the context of NR is the fact that the first three of 
these phases should be well described by classical and semi-classical calculations 
provided the BH mass exceeds the Planck scale by at least a factor of a fewP^D 
If we further assume that most of the energy of the collision process resides in the 
kinetic energy of the particles, such that their internal structure becomes negligible, 
the dynamics of the collision should be well modeled by two point particles or BHs 
in d-dimensional general relativity ESEHl Testing this assumption forms one of the 
main motivations for NR applications and we shall return to this question shortly. 

BH formation is expected to manifest itself in collision experiments by a special 
signature in its decay products as for example the jet multiplicity or transverse 
energyfHS F or the identification of these signatures, theoretical predictions from 
Monte-Carlo event generators such as BlackMaxPH and Charybdi^ 1 ^ 11 ^ 1 are 
compared with experimental data. Key input parameters for the event generators 
are the scattering cross section for BH formation and the initial mass and spin 
distributions of the formed holes. Providing this information forms a second main 
challenge for NR to complement estimates obtained from perturbative calculations 
on superposed Aichelburg-Sexl shock-wave backgrounds .UiZtfUH 
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The most recent analysis of data taken from collisions at the LHC with energies 
up to 8 TeV excludes semi-classical BHs of mass below 4.1 — 6.1 TeV^^ and sets 
exclusion contours in the plane formed by the BH mass and the threshold mass 
separating the semi-classical from the quantum regime It thus appears unlikely 
that the most optimistic TeV-gravity models with Mpi^ « 1 TeV are those realized 
in nature. 

In summary, the main questions to be addressed by NR simulations are (i) the 
validity of the Hoop conjecture in highly dynamical configurations, (ii) the scattering 
threshold for BH formation, (iii) the mass and spin parameters of the formed BHs 
and (iv) the impact of the colliding objects' structure on the dynamics. 

High-energy collisions in the framework of general relativity are currently best 
understood in d — 4 dimensions, largely because the numerical infrastructure for 
this setting is most advanced and robust. Even though generalization of the results 
to arbitrary spacetime dimensionality will ultimately be essential, a great deal of 
insight can be and has already been obtained in the four-dimensional case. We shall 
review these results before plunging into the latest developments on collisions for 
arbitrary d. 

The expectation of BH formation in high-energy collisions dominated by gravita- 
tional interaction is closely tied to Thome's Hoop conjecture.^ Applied to head-on 
collisions with boost factor 7, we expect two objects, each of rest mass mo and size 
i?o , to form a BH if the Schwarzschild radius R — 2 M = 477710 associated with M 
satisfies R > i? , i- e. 7 > i?o/(4m ). Choptuik & PretoriuJ^ tested this scenario 
by numerically simulating head-on collisions of boosted boson stars in d — 4 di- 
mensions. Their initial configurations consist of two boson stars with compactness 
2m n /R n w 1/20 and the Hoop conjecture would imply BH formation for boost fac- 
tors 7 > 10. Their numerical results demonstrate BH formation for 7 > 2.9 ± 10 %, 
a factor a few below the Hoop conjecture and thus strengthening the expectation 
that kinetic-energy dominated collisions indeed produce BHs. This picture has re- 
cently been confirmed in head-on collisions of fluid particles in Refs. 11571 11581 The 
energy radiated in gravitational waves (GWs) reported by East & Pretorius^^ is 
~ 16 ± 2 % for 7 = 10, in good agreement with ultrarelativistic limits 16.4 % and 
14 ± 3 % obtained from perturbative calculations^^"^! an d numerical simulations 
of BH head-on collisions'^ respectively. These studies strengthen the conjecture 
that trans-Planckian collisions are well described by high-energy collisions of black 
holes, general relativity's closest analog of point particles. 

Binary configurations of non-spinning, equal-mass BH binaries are characterized 
by two parameters, the boost factor 7 = — v 2 and the impact parameter 

defined as b = L/P, where v, L and P are the velocity, orbital angular momentum 
and linear momentum. The head-on case b = has been studied by Sperhake et 
al^SSi wn0 evolve a sequence of collisions with 7 < 2.93. The simulations always 
result in a merger and extrapolation of the radiated energy to v — > 1 leads to the 
above mentioned i? ra d = 14 ± 3 % of the Center-of-Mass (CoM) energy, about a 
factor 2 below Penrose's upper limit 1 — 1/V% quoted in ll61l . For grazing collisions, 
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the outcome depends on the size of the impact parameter b. Sperhake et a/P^I 
identify three regimes: (i) For b > b Bcat , the BHs scatter off each other and eventually 
escape to infinity, (ii) In a narrow regime b* < b < 6 sca t the BHs separate after an 
initial encounter, but have lost enough energy in gravitational radiation such that 
they form a bound system and eventually merge. As b approaches the threshold of 
immediate merger b* , the binary exhibits the zoom-whirl behaviour identified by 
Pretorius & KhuranaP^l w ith a number of whirls proportional to the logarithmic 
distance of b from b* . (iii) For sufficiently small b < b* , the outcome is a prompt 
merger. By fitting numerical results obtained for binaries with 7 < 2.3, Shibata et 
aiP^ have obtained a functional relation for the scattering threshold 

2. 5 ± 0.05, r 

Oscat = M 

V 

The numerical studies further demonstrate that grazing collisions with impact pa- 
rameters near the range b* < & scat generate enormous amounts of gravitational 
radiation up to at least ~ 35 % of the CoM energy; for a detailed analysis of the 
GW emission see Berti et a/PE3 Indeed it has been conjectured by Pretorius & Khu- 
rangP^H that all kinetic energy may be radiated in the form of GWs. This conjecture 
has recently been investigated numerically in Sperhake et aZP^ by monitoring the 
BH horizon properties in collisions of spinning and non-spinning binaries. While 
most of the kinetic energy is indeed radiated for mild boosts 7 < 1.5, GW absorp- 
tion by the holes becomes increasingly efficient for larger collision velocities; cf. their 
Fig. 1. Extrapolation of the numerical results to the ultrarelativistic limit predicts 
that approximately equal fractions of the kinetic energy are absorbed by the BHs as 
are lost in gravitational radiation, and sets a lower bound on the final BH mass of 
about half of the CoM energy. This study further demonstrates that the impact of 
the BH spin on the scattering threshold and radiated energy becomes negligible for 
7 > 2.5 lending additional support to the modeling of trans-Planckian scattering in 
terms of high-energy collisions of BHs. 

The most likely scenarios for BH formation at the LHC are proton-proton col- 
lisions. It will be important, therefore, to investigate the impact of electric charge 
on the dynamics which has been ignored in the above mentioned studies. A first 
exploration of colliding BHs with electric charge has been presented in Zilhao et 
a ^ji67j They extract gravitational as well as electromagnetic radiation generated in 
head-on collisions of BHs with equal mass and charge starting from rest. As in- 
tuitively expected, the electric repulsion slows down the dynamics and leads to a 
decreasing amount of GW energy Ef* d as the mass-to-charge ratio Q/M increases 
from towards the extremal case Q = M which corresponds to a static, albeit un- 
stable, configuration. The radiated electromagnetic energy Ef ad , on the other hand, 
is maximal near Q/M = 0.6 but always remains below the gravitational one with 
the quotient Ef^ d /E^ d increasing monotonically as a function of Q/M reaching 
0.25 as Q/M -> 1. 

In contrast to the four-dimensional case, BH collisions in d > 5 spacetime di- 



(7) 
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mensions appear to be more challenging from the view point of numerical stability 
and have not yet been explored to the same extent. Collisions of two equal-mass, 
non-spinning BHs starting from rest have been studied by Witek et a ;!i£MHID who 
use the Kodama-Ishibashi formalisrrP^ for the extraction of GWs and find that 
0.089 % of the CoM energy is radiated in GWs for d = 5, to be compared with 
0.05 % in four dimensions P^l These results were generalized to the case of unequal 
mass ratios r\ := mim-i/{mx + m-i) 2 < 0.25 in 11691 (cf. their Fig. 3) leading to the 
cxtreme-mass-ratio limit Ef* d /(Mr] 2 ) = 0.0164 - 0.0336ry 2 . 

Point-particle calculations are in good qualitative and quantitative agreement 
with these NR results^S but additionally make very interesting predictions for BH 
head-on collisions with non-zero initial velocitiesP^£^ d typically decreases with d 
to such an extent that for d = 11, head-on collisions starting from rest may radiate 
more GW energy than their ultrarelativistic counterparts; cf. Fig. 1 in 11711 This 
counter-intuitive prediction has yet to be tested with NR methods. BH collisions in 
d = 5 dimensional cylindrical spacetimes have been simulated in !172l The cylindrical 
setup can be viewed in terms of an infinite array of BHs resulting in a breaking effect 
on the collision dynamics. 

The most advanced study of BH scattering in d — 5 has been performed by 
Okawa et a/P^ who find scattering thresholds 6 sca t decreasing in terms of the 
Schwarzschild radius ro, = ^t 1 /^ -3 ) from 3.6 rs at v — 0.4 to 3.3 ro, at v = 0.6. 
Above this velocity, numerical instabilities prevent a determination of 6 S cat- In con- 
trast to the four-dimensional case, however, they do not observe any signs of zoom- 
whirl behaviour, probably due the equal fall-off ~ 1/r 2 of the gravitational and cen- 
trifugal potential in d — 5. By analysing the Kretschmann scalar C := R a ^ lS R a /3js , 
they further conclude that super-Planckian physics may be visible in the scattering 
of BHs in five spacetime dimensions. A straightforward calculation shows that for 
a single BH of mass M, the Kretschmann scalar on the horizon is given bj0 

By calculating \fC in the orbital plane, Okawa et al. observe super-Planckian do- 
mains where C S> Chor outside the individual holes horizons even for non-merging 
configurations. Regions with curvature radii below the Planck scale may thus form 
in high-energy particle collisions without being cloaked inside an event horizon. 

5. Gauge-gravity duality 

The Gauge-gravity duality conjectures the mathematical equivalence or duality be- 
tween field theories involving gravity in d dimensions on the one side and field 
theories without gravity in d — 1 dimensions on the other. It is one of the most 



b Note the difference of a factor 37r/8 arising from the different choices of G in the Einstein 
equations, our Eq. J2]l and Eq. (1) in !173l 
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remarkable results of string theory and, because of Maldacena's prototypical exam- 
p kj m H H6] of the duality between type lib string theory in 5 dimensional AdS space 
times the S 5 sphere and the J\f = 4 Supersymmetric SU(N) Yang- Mills (SYM) 
CFT in 4 dimensions, often referred to as the AdS/CFT correspondence. 

Following Maldacena's original example, a large number of similar dualities have 
been conjectured; see for example chapter 13 in Ref. 11771 Many of these involve type 
lib string theory or M theory compactified on AdS„ times a sphere or torus. The 
best understood example to date, however, is given by Maldacena's prototype which 
is also the version underlying present NR studies and, hence, the one we shall focus 
on in the remainder of this section. 

The AdS/CFT correspondence, while not rigorously proven, is inspired and sup- 
ported by a variety of features displayed by the two theories involved: 

• Holography: Black holes can be regarded as thermodynamic systems with 
entropy S = Ah OT /4G, where Ah m is the horizon area! 178 ^ 179 ! In contrast, 
the entropy of systems described by local field theories grows in proportion 
to the volume, suggesting a holographic connection to gravity EHflhiSlJ 

• Symmetries: The isometry group of AdS 5 x S 5 is 50(2, 4) x SO (6) which also 
leave Af = 4 SYM invariant. Fermionic symmetries agree similarly and lead 
to a supergroup SU(2, 2|4) which is identical to the JV = 4 superconformal 
symmetry; cf. Ref. 11831 

• D3 branes: A stack of N D2> branes has an open-string description equiva- 
lent to JV = 4 SYMp^ The stack also forms a solution of type lib super- 
gravity, i. e. the low-energy limit of type lib string theory^ 8 -^ which in the 
near-horizon limit of the extremal solution reduces to the AdS§ x S 5 met- 
ric. This equivalence further relates the coupling constant g and number of 
colors N of the gauge theory with the string length and coupling £ s , g s and 
the spacetime curvature radius L on the gravity side. Specifically, 

. 9 2 =47r ffs , 9 2 N=^. (9) 

• The Hawking-Bekenstein entropy associated with the AdS$ BH of the near- 
extremal Di brane is equal to the entropy of a non-ideal gas in JV = 4 SYM 
in the strong 't Hooft coupling limiP^ g 2 N -> oo; cf. [T53 H551 IT551 and 
Sec. 2.3 in[HS 

• The ratio of shear viscosity 77 to entropy density s is conjectured to have 
a lower bound given by the value rj/s = h/(An) which is universal to all 
theories with a gravity dual. 190 This bound appears to be satisfied by all 
known liquids. 

• AdS/CFT based studies (cf. below) predict rapid thermalization on time 
scales ~ 1 fm/c for quark-gluon plasma produced at Brookhaven's Rela- 
tivists Heavy Ion Collider (RHIC) in agreement with observations 



One of the most valuable properties of the correspondence arises from Eq. (j^J) . The 
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duality is strongly /weakly coupled and can be applied in both directions: strong 
coupling of the gauge theory g 2 N 3> 1 corresponds to L ^> £„, i. e. the regime of va- 
lidity of general relativity, while the gauge theory is rather well understood for small 
g 2 N and may provide information about the strongly coupled limit of string-theory. 
We need to bear in mind, however, that the correspondence applies to cousins of 
quantum chromodynamics (QCD) instead of QCD itself. Nevertheless, the corre- 
spondence provides an unprecedented opportunity to theoretically model physical 
systems in the strong-coupling regime of field theories and, as demonstrated by 
the above mentioned studies on shear viscosity and thermalization of quark-gluon 
plasma, is capable of making qualitatively and even quantitatively accurate predic- 
tions for physical systems. Reviews of the AdS/CFT correspondence with varying 



levels of detail are given in [1921 [1931 EEa MM H23 LLZ3 For an overview of poten- 
tial applications in fluid dynamics, condensed matter physics, super conductors and 
super fluidity we refer to Refs. QJH [Ml [Ml [Ml HH I20T1 

Anti-de Sitter space in d dimensions AdS^ is the maximally symmetric solution 
of the d-dimensional vacuum Einstein equations with cosmological constant A < 



R u 



A 9fW = 8irGT^ = . 



(10) 



It can be represented as the hyperboloid X$ - 
dimensional flat spacetime of signature 



Ei=i x ? embedded in a d + 1- 
with metric 



d-l 



ds A 



Transforming to coordinates 



-dX 2 



dXj 



(11) 



i=l 



COS T sin T 

X a = L , X d = L , Xi=Lt<mpni for i = l...d-l, (12) 

cos p cos p 

with hyperspherical coordinate^ Ylt=i = 1 , we obtain AdS f j in global coordi- 
nates with metric 

L 2 



ds 2 



, . -dr z + dp 2 - + s'm z pdfli 9 
cos 2 p K ' 



(13) 



where < p < n/2, —tt < t < tt. By unwrapping the cylindrical direction, the 
range of the time coordinate is often extended to r € R. 
Alternatively, Poincare coordinates (t,z,x l ) defined by 

+ L 3 + X>') 2 -! 2 



2z 



2z 



i=l 
d-2 



X, = , for i = l.,.d-2, 

z 



-L 2 + J2(x 1 ) 2 =t 2 



■X-d — ) 

z 



(14) 



c e. g. in d 
HQ 2 — 



= 5: Oi = sin x sin 6 cos </>, Q2 = sin x sin 9 sin <f>, Q3 = sin x cos 9, Q4 = cos\ and 
d X 2 + sin 2 X (d9 2 + sin 2 9d<t> 2 ). 
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give the AdS metric in the form 



T 2 



-dt 2 + dz 2 + ^2(d. 



d-2 

i\2 



X 

i=l 



(15) 



with z > 0, t E R. Poincare coordinates cover only half the hyperboloid (the other 
half corresponding to z < 0) and this patch is often referred to as the Poincare 
wedge. For its relation to global AdS see for example 202, 203, 204. 

At the AdS boundary p —> tt/2 or z — >• 0, the AdS metric becomes singular. 
It therefore induces a conformal class of metrics at the boundary, i. e. the induced 
metric is determined up to a conformal rescaling only. This remaining freedom 
manifests itself in the boundary topology of the global and Poincare metrics. In the 
limit p— >-7r/2orz— !>0, Eqs. (fT3|) and (|15|) become respectively 

d-2 

ds 2 gl ~ ~dr 2 + dQ 2 d _ 2 , ds 2 P ~ -dt 2 + d(x 1 ) 2 , (16) 

i=l 

which can be related by a conformal transformation; cf. Eq. (8.17) in 1 1771 In conse- 
quence, gravity in global or Poincare AdS is related to field theories on spacetimes 
of different topology: R x Sd-2 in the former and R d_1 in the latter case. 

There remains the extraction of physical quantities of the field theory from the 
gravitational side. This dictionary between bulk and boundary physics is estab- 
lished by the equivalence of the gravitational action of the bulk taken in the limit 
of the boundary and the effective action of the field theory on the boundary.E^E^ 
This duality involves divergences on both sides, infrared divergences on the gravity 
side due to the singular nature of the metric on the boundary and ultraviolet diver- 
gences of the quantum field theory. These anomalies and their elimination through 
renormalization have been discussed in 12051 12061 12071 NR applications have so far 
focused on the calculation of the vacuum expectation values of the field theory's 
energy momentum tensor (TYj) corresponding to the gravitational sector in the bulk 
and we shall only consider this case here. For a discussion including matter fields 
see for example Sec. 5 in 12061 as well as [207. 

Through the AdS/CFT correspondence, the vacuum expectation values (Ty) 
of the field theory are given by the quasi-local Brown- York^CSI stress-energy tensor 
and thus directly related to the bulk metric. Following de Haro et aZ.p^it is con- 
venient to consider the (asymptotically AdS) bulk metric in Fefferman-GraharrP^I 
coordinates 

L 2 

ds 2 — g^ u dx^dx v = — [dr 2 + ^fijdx 1 dx^\ , (17) 

where 

lij = 7ii(n x l = 7(0)»j + »~ 2 7(2)y + • • • + r d -i {d)l3 h {d)lJ r d logr 2 + 0{r d+1 ) . (18) 

Here, the "f( a )ij an d ^(d)ij are functions of the boundary coordinates x l , the loga- 
rithmic term only appears for even d and powers of r are exclusively even up to 



January 18, 2013 2:35 WSPC/INSTRUCTION FILE main 



16 Ulrich Sperhake 

order d—l. As shown in 206] the vacuum expectation value of the CFT momentum 
tensor for d = 4 dimensions is then obtained from 

^ ; 2 ..fem.jTi 
7(2) - 7(0) 7(0)7(2)fc;7(2)mn 

- \l{2)i m l{2)jm + J7(2)y7(2) I , (19) 

and 7(2)y is determined in terms of 7(o)ij- The dynamic freedom of the CFT is thus 
encapsulated in the fourth-order term 7(4)^. Note that for r — > 0, the metric ([T7|) 
asymptotes to the AdS metric in Poincare coordinates ([T5j) . The relation between 
the gravitational and CFT energy-momentum tensor for the case of global AdS is 
discussed in I205> see also 12101 

Many applications of the AdS / CFT correspondence are concerned with the equi- 
libration of matter in heavy-ion collisions at the RHIC or LHC and in particular its 
rapid thermalizationP^M^ While the quark-gluon plasma generated in the colli- 
sions is far from equilibrium at early stages, its behavior appears to be well described 
by hydrodynamics after time scales of the order 1 fm/c. This process is in principle 
governed by QCD but results indicate that many physical aspects can be studied 
in the framework of Af = 4 SYM through the gauge-gravity duality. Small pertur- 
bations of a static system in thermal equilibrium are known to decay exponentially 
fast and correspond to quasi-normal modes on the gravity sid e ! 214 ! and, as we shall 
discuss, numerical studies in the context of the AdS/CFT correspondence yield a 
similar picture for far-from equilibrium configurations. 

The use of NR methods to study these processes in AdS/CFT has been pio- 
neered by Chesler & Yaffe^J! 4 ^ who evolve an anisotropic source on the AdS bound- 
ary switched on after a short time using a characteristic approach based on ingoing 
Eddington-Finkelstein coordinates on the Poincare patch of AdS. Their numerical 
scheme is reduced to an effective "1+1" scheme by assuming boost invariance as 
well as rotational and translation symmetry in the transverse plane. Their boundary 
data generate gravitational waves which propagate into the bulk and lead to forma- 
tion of a BH. They extract the energy momentum tensor of the CFT dual and follow 
the time evolution of the energy density as well as the transverse and longitudinal 
pressure components. Isotropization of the pressure occurs on time scales inversely 
proportional to the local temperature at the onset of the hydrodynamic regime 
which translates into an isotropization time of 0.5 fm/c assuming a temperature 
of 350 MeV. Using the same setup, Chesler & Teaney^3 use different definitions 
of a "temperature" based on the energy density and on two-point functions. These 
two versions start agreeing after a time 1 fm/c coincident with the isotropization 
of transverse and longitudinal pressure. Wu & Romatschke^^ employ a similar ap- 
proach to collide two superposed Shockwaves in a boost invariant approximation 
and find that the late-time behavior of the energy density is given by a hydrody- 
namic description involving a scale parameter determined by the initial apparent 
horizon area. 
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A comparison of the fully non-linear numerical results with predictions from the 
linearized close- limit approximation^^ was performed by Heller et alW^ Instead of 
sourcing the anisotropy through a boundary term, they prescribe anisotropic data 
on an initial null hypersurface extending through the bulk. Their results confirm 
the short isotropization times ~ 1 fm/c and find the linear approach to reproduce 
these values within 20 % even for large initial anisotropies. 

A numerical scheme based on the ADM formalism of the Einstein equations was 
developed by Heller el a /!212J2lSl wn0 emphasize that thermalization, when defined 
as the onset of the hydrodynamical regime^, may differ from isotropization. By 
evolving boost-invariant, transversely homogeneous plasmas with various different 
initial conditions, they find that a hydrodynamic description may well be applicable 
when the pressure is still anisotropic. 

In Ref. |43j Chesler & Yaffe relax their symmetry assumptions to translation 
invariance in the transverse direction which effectively constitutes a 2+1 numerical 
scheme. Their characteristic formulation works well in this scenario without showing 
any signs of formation of caustics and enables them to model heavy-ion collisions by 
colliding two shock waves. Specifically, they consider a single shock-wave solution in 
Fefferman- Graham coordinates, superpose two of those and transform the result to 
ingoing Eddington-Finkelstein coordinates. For stability purposes, they introduce 
a small energy offset on the CFT side which generates an apparent horizon in the 
gravity dual above the Poincare boundary to absorb steep gradients encountered in 
the metric functions deep in the bulk. By comparing the numerically determined 
pressure components with hydrodynamical predictions, they confirm the picture of 
rapid isotropization obtained in scenarios of higher symmetry. Translated into values 
for gold ion collisions at the RHIC, they observe isotropization about 0.35 fm/c after 
their shock waves start overlapping. 

A numerical code based on Cauchy type evolutions of asymptotically global AdS 
spacetimes using the GHG scheme has been developed by Bantilan et alW^ By as- 
suming an 5*0(3) symmetry, they reduce their computational domain to 2-1-1 dimen- 
sions. Initial data is specified in the form of a localized scalar field which promptly 
collapses to a BH with a highly distorted horizon and settles down to a stationary 
configuration through quasi-normal ringdown. Whereas the lowest ringdown modes 
agree with linearized predictions, higher angular modes exhibit significant coupling 
due to non-linear effects. The dual stress-energy tensor of the CFT is mapped from 
the global RxS 3 AdS boundary to a Minkowski background and found to evolve in 
agreement with that of a thermalized SYM fluid from the start of the simulations. 
Their work furthermore discusses in detail a number of regularization procedures 
required to obtain a numerically stable framework. 

A very recent paper by Adams et a/P22 presents a first numerical investigation 
of superfluid turbulence in 2+1 dimensions corresponding to a 3+1 dimensional 
asymptotically AdS gravity dual. They observe energy injected at long wavelengths 



d see Eq. (5) in Ref. 12191 for their precise definition 
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to cascade down to short length scales, in contrast to the behavior of normal fluid 
turbulence in 2+1 dimensions. 

6. Conclusions 

The vast majority of the work reviewed here has been generated during the last four 
years; NR in higher dimensions is a very young research field and likely has only 
scratched the surface of a wealth of possible applications. Nevertheless, a number 
of impressive results have demonstrated its vast potential. Discoveries such as the 
formation of naked singularities in black strings strengthen the belief that higher- 
dimensional gravity differs from its four-dimensional counterpart in more than a 
quantitative manner. The same holds for the asymptotic structure of the space- 
times under consideration. The instability of asymptotically anti-de Sitter space- 
times to BH formation under arbitrarily small perturbations is in stark contrast to 
the stability of asymptotically flat or de Sitter spacetimes 1222023 

One of the remarkable properties of higher-dimensional gravity is its connection 
with other areas of physics. The holographic principle and the AdS/CFT correspon- 
dence in particular have opened a path towards understanding a variety of high- 
energy physics phenomena through the study of 4+1 gravity. This facilitates direct 
contact with experiments such as heavy- ion collisions performed at the RHIC. At 
the same time, parton-parton collisions at TeV energies at the LHC have provided 
for the first time a direct experimental test of the existence of extra dimensions. 

In spite of the great success of NR in recent years, a great deal of work remains to 
be done. Numerical studies in d dimensions have as yet been restricted to systems of 
considerable symmetry. This is largely a consequence of computational requirements 
which increase enormously with each extra dimension. Continuous improvements of 
computer technology should, however, make possible in the foreseeable future the 
study of generic systems at least in d — 5, opening up a much wider class of problems 
in the context of the AdS/CFT correspondence and the stability analysis of BHs. 

Even within the context of symmetric systems, many challenges are awaiting the 
community. Trans-Planckian scattering is rather well-understood in d = 4. Connec- 
tion with experiment, however, requires extension of these results to higher d. Here, 
progress has been obstructed by stability issues of the numerical codes, calling for 
further investigation of the numerical formulations and gauge conditions. Many 
open questions also remain about the stability of BHs in higher dimensions. Here 
we note in particular that BHs are no longer restricted to spherical topology and 
even stationary multiple-BH solutions are known. It is unclear how present numer- 
ical frameworks are able to handle black rings or saturns, although the relatively 
straightforward extension of NR studies to black strings may be encouraging. Nu- 
merical codes developed for the AdS / CFT correspondence have as yet been tailored 
to rather specific problems. It would be desirable to have a generic framework ef- 
fectively representing a numerical AdS/CFT laboratory. Given the enormous range 
of BH systems in the context of the correspondence, this is likely going to be a 
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tremendous challenge, albeit one of large scientific potential. 

In summary, numerical relativity in higher dimensions is a highly active and 
rapidly growing held and it will be very exciting to see its development over the 
next years. 
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